f{R) term and geometric origin of the dark sector in Horava-Lifshitz gravity 
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Inclusion of ,f{R) term in the action of Horava-Lifshitz quantum gravity with projectability but 
without detailed balance condition is investigated, where R denotes the 3-spatial dimensional Ricci 
scalar. Conditions for the spin-0 graviton to be free of ghosts and instability are studied. The 
requirement that the theory reduce to general relativity in the IR makes the scalar mode unstable 
in the Minkowski background but stable in the de Sitter. It is remarkable that the dark sector, dark 
matter and dark energy, of the universe has a naturally geometric origin in such a setup. Bouncing 
universes can also be constructed. Scalar perturbations in the FRW backgrounds with non-zero 
curvature are presented. 
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I. INTRODUCTION 

There has been considerable interest recently on a the- 
ory of quantum gravity proposed by Horava [l|, moti- 
vated by Lifshitz scalar theory 0. A free scalar field in 
(d + l)-dimensional flat spacetime is usually given by 



dt(fx 



(1.1) 



which is Lorentz invariant, 

i->e°(t,x), x'^e(t,x), (z,j = 1,2,3). (1.2) 

The corresponding propagator is given by 6*0(0;, k) — 
1/fc^, indicating that the theory is not renormalizable. 
To improve the ultra-violet (UV) behavior, Lifshitz in- 
troduced high-order spatial derivatives, (j) (— V^) </>, into 
the action, and found that the resulted theory becomes 
renormalizable for z > d [3]. An immediate consequence 
of these terms is that the theory is no longer Lorentz 
invariance, and t and x scale differently, 



t Pt, x' ix' 



(1.3) 



Based on the above observations, Horava toke the point 
of view that Lorentz symmetry should appear as an emer- 
gent symmetry at long distances, but can be fundamen- 
tally absent at short ones 4, 5]. To realize such a per- 
spective, Horava started with the Arnowitt-Deser-Misner 
(ADM) form of the metric. 



ds' 



N^c^dt^ + g,j {dx' + N'dt) {dx^ + N^dt) , 



(1.4) 



and imposed the foliation-preserving diffeomorphisms. 



(1.5) 



to be denoted by Diffjr(M). At low energies, the theory 
is expected to flow to the IR fixed point z — I, whereby 
the Lorentz invariance is "accidentally restored." 

The effective speed of light in this theory diverges in 
the UV regime, which could potentially resolve the hori- 
zon problem without invoking inflation Q. The spatial 



curvature is enhanced by higher-order curvature terms 
0-0 , and this opens a new approach to the flatness prob- 
lem and to a bouncing universe 

BiEini- In addi- 
tion, in the super-horizon region scale- invariant curvature 
perturbations can be produced without inflation fl^-^T?] . 
The perturbations become adiabatic during slow-roll in- 
flation driven by a single scalar field, and the comoving 
curvature perturbation is constant [l7[. Due to all these 
remarkable features, the HL theory has attracted lot of 
attention lately, see, for example, Ref. [T^ and references 
therein. 

To formulate his theory, Horava started with two con- 
ditions - detailed balance and projectability [l[. The 
detailed balance condition restricts the form of a gen- 
eral potential in a (d + l)-dimensional Lorentz action 
to a specific form that can be expressed in terms of a 
d-dimensional action of a relativistic theory with Eu- 
clidean signature, whereby the number of independent- 
couplings is considerably limited. The projectability con- 
dition, on the other hand, restricts the lapse function N 
to be space-independent, while the shift vector N"^ and 
the 3-dimensional metric gij still depend on both time 
and space. 



N = N{t), = N\t,x), - .g„-(t,x). 



(1.6) 



Clearly, this condition is preserved by Diffjr(M). 

It should be noted that, due to the restricted Diffjr(M), 
one more degree of freedom appears in the gravitational 
sector - a spin-0 graviton. In particular, in the pro- 
jectable case this mode is unstable in the Minkowski 
background [1, [l^, [2l| . As shown below, this is a generic 
feature of the theory with projectability condition. How- 
ever, this instability does not necessarily show up in phys- 
ical environments |22l |. 

In addition, it is also possible that the legitimate back- 
ground in the HL theory is not Minkowski ^. In par- 
ticular, recent observations show that our universe is 
currently de Sitter-like Therefore, instead of the 



^ Thanks to P. Horava for pointing it out to us. 
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Minkowski, one may take the de Sitter space as the back- 
ground. As a matter of fact, it was shown recently that 
the de Sitter space is indeed stable in the framework of 
the HL theory with projectability condition [l8l[T9|. 

On the other hand, in the non-projectable case, this 
mode is also unstable [2J]- However, if one includes terms 
made of the spatial gradients of the lapse function. 



example, in [33| TZ was token as 



^ainiV, 



(1.7) 



the mode can be rendered stable [2g] . However, an imme- 
diate price to pay for this is the inclusion of more than 60 
sixth-order derivative terms in the potential [2^. When 
matter, such as a scalar field, is included, the number 
of such terms increases dramatically. In addition, strong 
coupling may still exist 27[, unless the scales appearing 



in hi gher order terms are much lower than the Planck 
scale |25| . 

It should be noted that strong couplings also happen 
in the version with projectability condition [28j . How- 
ever, Mukohyama recently showed that, when nonlinear 
effects are taken into account, the spin-0 graviton decou- 
ples for spherically symmetric, static, vacuum spacetimes 
|29| . Together with Wu, we showed that this is also the 
case in the cosmological setting [19]. As a result, the rel- 
ativistic continuation A — 1+ exists in the IR, whereby 
the strong coupling problem is circumvented. 

The most general form of the HL theory with pro- 
jectability but without detailed balance condition was 
first developed in [13, HH by Sotiriou, Visser and We- 
infurtner (SVW), in which the highest order of spatial 
derivatives is assumed to be six, the minimal require- 
ment for the theory to be power-counting renormalizable 
in (3-1-1 )-dimensional spacetimes [Ull]- However, in prin- 
ciple there is nothing to prevent one to construct actions 
with spatial derivatives higher than six. In addition, in 
condensed matter physics, the critical exponent z is also 
not necessarily an integer. In this paper, we explore all 
these possibilities. Clearly, promoting all terms to high 
orders is out of control. Instead, we consider high or- 
der derivatives only coming from the Ricci scalar R, by 
simply replying the third-order polynomial of R in the 
SVW setup by an arbitrary function, f{R), while keep 
the rest the same. As shown below, this simple gener- 
alization results in very rich physical phenomena, and 
in particular can give rise to dark matter and dark en- 
ergy. Since in such a setup the origin of them is purely 
geometric, it automatically explains why these objects 
are "dark" . Bouncing universes can be also easily con- 
structed by properly choosing the form of f{R). 

It is interesting to note that in [30| it was advocated 
that the HL theory with projectability condition has a 
built-in dark matter component, due to the non- locality 
of the Hamiltonian constraint. In addition, f{TZ) models 
have been investigated in the framework of the HL theory 
[sH - issj . but in all these studies TZ is different from the 
3D Ricci scalar R. As a result, those f{TZ) models are 
fundamentally different from the ones studied here. For 



TZ = KijK'^ - \K'^ + 2^V^ [n^y^n" - n^V ^vP) 

~E'^g,,uE''\ (1.8) 

where (j/ = 0, 1, 2, 3) is a unit vector perpendicular to 
the hypersurfaces oit — Constant, A and [i are two cou- 
pling constants, Qiju is the"generahzed" De Witt metric, 
and E^^ is given in terms of a super-potential. For de- 
tail, we refer readers to (3ll - [35j . Clearly, for any choice 
of A, and i?'-' , due to the presence of the first term on 
the right-hand side of Eq. (|1.8p . TZ cannot reduce to the 
three-spatial dimensional Ricci scalar R. In addition, due 
to its presence, the corresponding theory usually involves 
high order time derivatives |3ll435j . a situation that was 
avoided in the first place by Horava, in order to circum- 
vent the ghost problem 1]. 

The rest of the paper is organized as follows: In Sec- 
tion II, we present the Hamiltonian and momentum con- 
straints, the dynamical equations, and the conservation 
laws of energy and momentum, after including the /(i?) 
term in the SVW setup. To separate the effect of this 
term from others introduced in [lO, [2l| , in this section 
we only consider the /(i?) term. In Section HI, we study 
the stability of the spin-0 graviton in the Minkowski back- 
ground, and obtain the stability condition. We show ex- 
plicitly that it is the requirement that the theory reduce 
to its relativistic limit at the IR leads to the instability 
of the spin-0 graviton. This is also true when all the 
other terms are included. In Section IV, we investigate 
its applications to cosmology, by first writing the corre- 
sponding Friedmann equations of the FRW universe with 
arbitrary curvature, and then their linear scalar pertur- 
bations. From these expressions, it can be easily shown 
that the de Sitter spacetime is stable, similar to the case 
without the /(i?) term [l^. In this section, we also gener- 
alize our studies to include all the other terms introduced 
in [20I . [2]| . In Section V, we present our main results and 
give some concluding remarks. 



II. INCLUSION OF THE f{R) TERM 

To see the role that the term /(i?) may play, we shall 
first neglect all the other terms constructed in j2^, i21j] , 
and simply consider the action, 



s = (- 



J dt(fxNy^{CK 



f{R) + C'C 



M 



where g = det < 



C 



K 



(2.1) 



(2.2) 



= I/IGttG, a is a dynamical coupling constant, and 
f{R) is an arbitrary function of R. However, to have the 
theory power-counting renormalizable, it is necessary to 
include terms that are equal or higher than R^. The 
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extrinsic curvature Ka is defined as 



K^, = ^ {-g^, + y^N, + V,iV,;) , (2.3) 

where Ni = gij , and the covariant derivatives refer to 
the three- metric gij. 

Variation with respect to the lapse function N{t) yields 
the global Hamiltonian constraint, 

J <fx^ {Ck - .f{R)) = SvrG j J* , (2.4) 



where 



J* = 2 



5{NCm) 
SN 



(2.5) 



Variation with respect to iV* yields the super-momentum 
constraint, 

VjTT^J' = 8ttGJ\ (2.6) 
where tt*-' and J* are defined as 

= -K'^ + \Kg'\ 



r = -N 



'SN~' 



(2.7) 



Varying with respect to g^ , on the other hand, leads to 
the dynamical equations. 



2XKK'^ ~2(K 



1 



+ — Vfc {N'^tt'^ - tt'^'N^ - n'^m') + SttGt'^ (2. 



where {K^Y^ = X^'if/, and 
^ 2 S{^Cm) 

F'l EE V'\7^F ~ g'^\/^F ~ FR'^ + ^fg'^, (2.9) 

with F = df{R)/dR. 

The matter quantities (J*, J*, r'-') satisfy the conser- 
vation laws, 



/ 



0, (2.10) 



^(V5J.)--^v.j'= 

~ (VfeTV, ~ V.ATfc) = 0. (2.11) 



III. SPIN-0 GRAVITON IN MINKOWSKI 
BACKGROUND 

It can be shown that Minkowski spacetime ds^ = 
—dt^ + Sijdx^dx^ is a solution of the vacuum field equa- 
tions, provided that /(O) ~ 0. Considering its linear 
scalar perturbations. 



5g.,j = -2V% -I- 2E^,j, SN, = B„ SN = 0, (2.1) 

in the quasi- longitudinal gauge (j) = Q = E, we find 
that to second order the action without matter takes the 
form. 



di^d^xN^g{c%^ +J^^\R)) 

- 3A) (3V? + 2^925 



/ dtd^x 
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+ (1 - \)Bd^B + - 8wC"V5V , (2.2) 

where 

7^-/'(0), ^^-CV"(0). (2.3) 
Variations with respect to B and ip yield, respectively, 
(1 - 3A)V'fc ~ (1 - \)k'Bk = 0, (2.4) 
(3A - 1) (i^k - k'Bk) = 2fc2 + Aojk'/Q') V'fe,(2.5) 



in the momentum space. When A 7^ 1, from Eq. 
we can express Bk in terms of tpk, and then Eq. 
becomes 



C(2) 



2e 



with 



,2 



dtd^x 



A-1 
3A - 1' 



c'fc' (7 + 



(v:;2_^2^2^, (2.6) 



4a; fc" 



,2 „ 



(2.7) 



Hence, Vfc satisfies i})k + w^V'fc = 0; which has stable 
solutions only when w^, > 0. Clearly, this is true in the 
IR limit only when 



A > 1, 7 > 0. 



(2.8) 



(2) 

These conditions also assure the kinetic term in Sg al- 
ways positive, that is, free of ghosts. When A = 1, from 
Ea. (|2.4p we find that ■i/'fc and Bk are independent of <, and 
so are the two gauge-invariant quantities ^ $fe — Bk and 
^fc = ■i/'fc. Therefore, the spin-0 graviton is stable when 
A = 1. 

Note that the addition of other terms presented in 
[20I [2H does not change the IR behavior, and their con- 
tributions only change the expression u'^ to the form, 



ctk' 7 + 



(4a;-353)fc' , (3^8 - 857) fc 



(2.9) 
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where gi are the couphng constants introduced in |2Cll21j . 
In addition, the condition that the theory reduces to gen- 
eral relativity (GR) requires 7 = — 1, a condition that 
was assumed in ^20, [2lj . This explains why in the S VW 
setup the spin-0 mode is not stable ^] . Our analysis given 
above shows that the instability is a generic feature of the 
HL theory with projectability condition. However, this 
instability does not necessarily show up in a physical en- 
vironment, as longer as some conditions are satisfied (2^ . 
In particular, if L/\c^\ > tj, the instability will not show 
up, where L is the length scale of interest, and tj denotes 
the timescale of Jeans instability [2^. In addition, the 
linear instability is stabilized by higher derivative terms 
if |c^| < l/(LAf*), where Af* is the energy scale sup- 
pressing higher derivative terms. In the current setup, it 
is given by ~ \g3\-^/^Mpi, \gs>5\-^/^Mpi. The lin- 
ear instability can be also tamed by Hubble friction, if 
L/|C^| > 1/H, where H is the Hubble expansion rate 
at the time of interest. Therefore, even a solution is not 
stable, one may properly choose the form of the function 
f{R), so that at least one of these conditions is satisfied. 
Then, the corresponding instability will not show up. 

It is also interesting to note that, although the SVW 
setup is not stable in the Minkowski background, it is 
stable in the de Sitter [31 • This is also true in the current 
setup. Then, one may consider the instability found in 
the Minkowski background is just an indication that the 
latter is not the ground state of the theory. Instead, the 
legitimate background in the HL theory might be the de 
Sitter spacetime. This point of view is further supported 
by current observations that our universe is currently de 
Sitter- like [H. 



The conservation of momentum (|2.11l) . on the other 
hand, is also satisfied identically, while the conservation 
of energy (|2.10|) gives 



'p + 3H (p + p) =0, 



(4.4) 



which can be obtained from Eqs. ()4.2|) and ()4.3|) . When 
K = we have R = 0, and these two equations reduce 
exactly to these given in GR with a cosmological constant 
A = -/(0)/2. In addition, Eqs. g^) and gS]) also 
admit a de Sitter point that corresponds to a vacuum 
solution {p = p = 0), at which /(i?) - Rf'{R) = -2A, 
which has the solution f{R) — — 2A — -fR. 

In the rest of this paper, we shall consider only the 
cases where k 7^ 0. Then, Eqs. ()4.2|) and (|4.3|) are 
quite different from the ones of 4-dimensional /(i?^"*^) 
models [s^ [s^l- In particular, since contains sec- 
ond derivatives of a, the generalized Friedmann equa- 
tions are usually high-order differential equations in the 
4-dimensional /(i?^^') models [sj- But, in the HL 
theory R[= 6k /a^] is a function of a only. As a result, 
/(i?) is a polynomial of a only, too. Therefore, in the 
current setup these terms act like sources. For example, 
if f{R) oc i?^/^, then from Eq. (|4.2|) we can see that this 
corresponds to dark matter. Note that to have F{R) 
real, in this case we must assume k > 0. In addition, for 
f{R) oc i?" Eq. g2|) shows that cx a^^". Then, with 
n < 1 this term can mimic dark energy. On the other 
hand, for n > 2 it can give rise to a bouncing universe 
0, H, 113, Ell • Since all of these terms are purely geomet- 
ric, they do not subjected to the energy conditions |3l| 
and matter instabilities. Note that the unification of the 
dark sector in /(^''^i?) models was first considered in [39|. 



IV. COSMOLOGICAL MODELS 



B. Linear Scalar Perturbations 



A. FRW Background 

The homogeneous and isotropic universe is described 
by the FRW metric, ds^ — ~dt^ + a?{t)^ijdx^dx^ where 

7ij = (1 -|- ^Kr^) with k — 0, ±1. For this metric, 

Kij = —a?H^ij and Rij — 2^7^^, where H — aja and an 
overbar denotes a background quantity. Setting Q 



-2p, r = 0, % 



(4.1) 



where p and p are the total density and pressure, we find 
that the Hamiltonian constraint (12.41) reduces to 



SttG 



(4.2) 



where R = 6k/ a^, and A = (3 A — 1) /2. It can be shown 
that the momentum constraint (|2.6p is satisfied identi- 
cally, while the dynamical equations (|2.8p yield, 



A^ 
a 



AttG 



(p + 3p) 



lfiR) + ^F{R). 



(4.3) 



Consider scalar perturbations in the quasi-longitudinal 
gauge [a]. 



SN = 0, SNi = a^Bu 



-2a^i/'7y: 



(4.5) 



it can be shown that the Hamiltonian constraint (12. 4p to 
first order yields, 

J y/^d^x F{R) + SkV') - AH (y^B + 3V/) 

= 47rGa2 J y/^d^x6p,, (4.6) 

where Sp. = —SJ*/2, V^/ = Y''f\ij: and \i denotes the 
covariant derivative with respect to "fij . The momentum 
constraint (12.61) takes the form. 



(3A ~l)ip' - 2kB + (A - 1) V^S = SnGaq, (4.7) 

with S,P = a^^gl*. The dynamical equations p.Sp . on 
the other hand, yield. 
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8F'(i?) 



3(3A- 1) 
iv^ (B' + 2HB) = ^SV, 



B' + 2-HB = i^(iT')V' 



AF'{R) 



STrGa^n, 



(4.8) 



(4.9) 



where 



{5V + 2p^) 7»J' +ni<'^'> 



ni<y> ^ ni'^' — 7'^v2n. 

3 ' 

The conservation laws give, 



q' + iUq - aiJT' - y (V^ + 3fc ) n = 0. 



(4.10) 

= 0, 
(4.11) 
(4.12) 



Following [18[ it can be shown from Eq. (|4.8p that the de 
Sitter background is indeed stable in the current setup. 



C. Inclusion of Other Terms Upto Six-Order 

As mentioned above, one can add the terms [1^ [2l| 



+ ^ [-fiRV^R + 75 (V,i?,fe) {V'W'')] , (4.13) 

into (|2.ip . so that up to the sixth-order derivatives the 
action is the most general one, where 7s's are the coupling 
constants, defined as {7s} = (53, ffs, 56j 57j Ss)- Then, the 
Hamiltonian constraint can be obtained from Eq. 
with the replacement of f{R) by — ^Cy, where 



Cv = 5Cv-f{R). 



(4.14) 



while the super-momentum constraint (j2.6p and the con- 
servation laws, (I2.10p and (|2.1ip remain the same. The 
dynamical equations are also given by Ea. (|2.8p . but now 
with 



(4.15) 



s=l 



where Fj|^^ is defined by Eq. jlH), ra^ = (1,2,2,2,2), 
and {J'sT' = {F^F5,Fe,Fr,F8y^, where {Fsf^'s are de- 



fined in d iS [23 

Due to the inclusion of the above terms in the action, 
Eqs. (|4?2l) and (|4?3| become 



AiJ^ = — p--^/(i?) + ^ + ^,(4.16) 



~a A'kG , , 1 „,^^^ K 

A- = ^ P + 3p - -/ i? + -F i? 

a 3 6 



where 



o _ 71 _ 372 73 

Pi = ^, = ^3 . 



(4.17) 



(4.18) 



Similarly, one can obtain the corresponding scalar per- 
turbations from the above and [9|] , which will not be pre- 
sented here. 



V. CONCLUSIONS 

We studied the efforts of spatial derivatives higher than 
six, and the possibility of the critical exponents z being 
non-integer, in the HL gravity with projectability but 
without detailed balance condition, by replying the third- 
order polynomial of the Ricci scalar in the SVW setup 
US [Hi by an arbitrary function, /(i?), while keeping the 
rest the same. The requirement that the theory reduce 
to GR in the IR makes scalar perturbations unstable in 
the Minkowski background, but stable in the de Sitter. 
This might imply that the legitimate background in the 
HL theory is not Minkowski spacetime, but the de Sitter. 

This simple generalization leads to very rich physical 
phenomena. In particular, it naturally gives rise to dark 
matter and dark energy. We also studied scalar perturba- 
tions in the FRW backgrounds with non-zero curvature, 
and presented the general formulas, which will make the 
studies of perturbations in such a setup considerably eas- 
ier. It would be very interesting to study their applica- 
tions to large-scale structure formation, early universe, 
and investigate constraints of the models from solar sys- 
tem tests. 

When high-order curvature terms, /(^''^i?), are in- 
cluded, it was shown that spin-2 massive models also 
exist and are ghosts [ioj . In the SVW setup of the HL 
theory, we have shown recently that vector perturbations 
vanish, while the only remaining tensor modes are the 
massless spin-2 modes [4l|. Therefore, it would be very 
interesting to see the effects of the f{R) terms on the 
vector and tensor modes in the current setup. 

In addition, dark matter and dark energy have been 
studied extensively [42| - [43 |. and current observational 
data impose serve constraints on various models [2^. It 
is also very interesting to fit the data to our current mod- 
els by using the MCMC code, recently developed by us 
[45I , including the studies of the growth factor of pertur- 
bations [46j . 

Strong couplings [11] may exist here, too. One way 
is to provoke the Vainshtein mechanism [47], similar to 
what was done in [l^ [2^ . 

Acknowledgements: The author would like to thank 
P. Horava, R. Maartens, A. Papazoglou, and D. Wands 



6 



for valuable suggestions and comments. Part of this 
work was done when the author attended the advanced 
workshop "Dark Energy and Fundamental Theory," sup- 
ported by the Special Fund for Theoretical Physics from 
the National Natural Science Foundation of China by 



the grant 10947203. The author would like to express his 
gratitude to Jianxin Lu for his organization of the Con- 
ference and hospitality. This work was supported in part 
by DOE Grant, DE-FG02-10ER41692. 



P. Horava, P hys. Rev. D79, 084008 (2009) 
[arXiv:0901.3775] . [26 
E.M. Lifshitz, "On the Theory of Second-order Phase 
Transitions I," Zh. Eksp. Toer. Fiz. 11, 255; 269 (1941). [27 
M. Visser, Phys. R ev. D80, 025011 (2009) 
!arXiv:0902.0590l ; larXiy:0912.4757 [28 
T.G. Pavlopoulos, Phys. Rev. 159, 1106 (1967). 
S. Chadha and H.B. Nielsen, Nucl. Phys. B217, 125 [29 
(1983). 

E. Kiritsis and G. Kofinas, Nucl. Phys. B821, 467 (2009) [30 
[arXiv:0904.1334 . 

G. Calcagni, J.' High Energy Phys., 09, 112 (2009) [31 
^Xiv:0904.0829 . 

G. Calcagni, Phys. Rev. D81, 044006 (2010) [32 

[arXi^:0905.3740| . 

A. Wa ng and R. Maart ens, Phys. Rev. D 81, 024009 [33 
(2010) arXiv:0907.1748;. 

R. Brandenberger, arXiv:0904.2835l [34 

A. Wang and Y. Wu, JCAP, 07, 012 (2009) 
[arXiv:0905.4117 . [35 
S. Mukohyama, J. Cosmol. Astropart. Phys., 06, 001 

(2009) arXiv:0904.2190 . [36 
S. Mukohyama, Phys. Rev. D80, 064005 (2009) 

farXiv:0905.3563 . [37 

Y.-S. Piao, Phys. Lett. B681, 1 (2009) larXiv:0904.4117l. 
R.-G. Cai, B. Hu, and H.-B. Zhang, Phys. Rev. D80, [38 
041501 (2009) arXiv :0905.0255 . 

B. Chen, S. Pi, and J.-Z. Tang, J. Cosmol. Astropart. [39 
Phys. 08, 007 (2009) ;arXiv:0905.2300 . 

A. Wang, D. Wands, and R. Maartens, J. Cosmol. As- [40 
tropart. Phys., 03, 013 (2010) arXiv:0909.5167 . 
Y.-Q. Huang, A. Wang , and Q. Wu, Mod. Phys. Lett. 
A25, 2267 (2010) [arXlV :1003.2003 . [41 
A. Wang and Q. Wu, "Stability of spin-O graviton and 
strong coupling in Horava-Lifshitz theory of gravity," [42 
arXiv: 1009.0268, 

T.P. Sotiriou, M. Visser, and S. Weinfurtner, Phys. Rev. [43 
Lett. 102, 251601 (2009) iarXiv:09 04.4464 . 

T.P. Sotiriou, M. Visser, and S. Weinfurtner, J. High [44 
Energy Phys., 10, 033 (2009) arXiv:0905.2798 . 
K. Izumi and S. Mukohyama, Phys. Rev. D 81, 044008 [45 

(2010) arXiv:0911.1814 . 

E. Komatsu, et al, arXiv:1001.4538 , [46 

C. Bogdanos, and E.N. Saridakis, Class. Quant. Grav. 
27, 075005 (2010) arXiv:0907 .1636] . [47 

[25] D.Blas, O.Pujolas, and S. Sibiryakov, Phys. Rev. Lett. 



07, 



197 



11, 



064028 



078 (2009) 
(2010) 



044004 (2010) 



185021 



104, 181302 (2010) arXi v:0909.3525] . 
L Kimpton and A. Padilla, J. High Energy Phys.. 
014 (2010) arXiv :1003.5666] . 

A. Papazogiou and T. P. Sotiriou, Phys. Lett. B685 
(2010) arXiv:0911.1299;. 

K. Koyama and F. Arroja, J. High Energy Phys., 03, 061 
(2010) arXiv:0910.1998 . 

S. Mukohyama, Class. Quant. Grav. 27, 223101 (2010) 
[arXiv: 1007.5199 . 

S. Mukohyama, j. Cosmol. Astropart. Phys., 0909, 005 

(2009) arXiv:0906.5069 . 
J. Kluson, J. High Energy Phys., 
arXiv:0907.3566 . 

J. Kluson, Phys. Rev. D81, 
arXiv:0910.5852 . ' 

J. Kluson, Phys. Rev. D82, 
arXiv: 1002.4859 . 

M. Chaichian et al, Class. Quantum Grav. 27, 

(2010) arXi v:1001.4102i . 

S. Carioni, et al, Phys. Rev. D82, 065020 (2010) 
[arXiv: 1003.3925 . 

T.P. Sotiriou, an d V. Fa raoni, Rev. Mod. Phys. 82, 451 
(2010) arXiv:080 5T726] . 

A. de Felice and S. Tsujikawa, Living Rev. Rel. 13, 3 
(2010) arXiv: 1002.4928 . 

S. Hawking and G.F.R. Ellis, The Large Scale Structure 
of Space- Time (Cambridge University Press, 1972). 
S. Nojiri and S. D. Odintsov, Phys. Rev. D68, 123512 
(2003) hep-th/0307288 . 

C. Bogdanos, S. Capozziello, M. de Laurentis, 
and S. Nesseris, Astropart. Phys. 34, 236 (2010) 
arXiv:0911.3094 . 

A. Wang, "Vector and tensor perturbations in Horava- 
Lifshitz cosmology," arXiv: 1008.3637 
A.M. Szelc, "Dark Matter Experimental Overview," 
arXiv:1010.3918 ; and references therein. 
S. Capozziello and M. Francaviglia, Gen. Relativ. Grav., 
40, 357 (2008) arXiv:0706 .1146j . 

S. Tsujikawa, "Drak Energy: investigation and model- 
ing," larXiv: 1004. 1493 and references therein. 
Y.-G. Gong, Q. Wu, and A. Wang, Astrophys. J. 681, 
27 (2008) arXiv:0708.1817 . 

Y. Gong, M. Ishak, A. Wang, Phys. Rev. D 80, 023002 
(2009) arXiv:0903.0001 . 

A. I. Vainshtein, Phys. Lett. B39, 393 (1972). 



